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_ Abstract 

^ . The classic Michaelis-Menten equation describes the catalytic activities for ensembles of enzyme 

C^ , molecules very well. But recent single-molecule experiment showed that the waiting time distribu- 
tion and other properties of single enzyme molecule are not consistent with the prediction based 

Ph ' on the viewpoint of ensemble. It has been contributed to the slow inner conformational changes 

o, 

, ^, . of single enzyme in the catalytic processes. In this work we study the general dynamics of single 

enzyme in the presence of dynamic disorder. We find that at two limiting cases, the slow reaction 

> 

Tij- ■ and nondiffusion limits, Michaelis-Menten equation exactly holds although the waiting time distri- 

\o _ 

en . bution has a multiexponential decay behaviors in the nondiffusion limit. Particularly, the classic 
'n|- ■ 

^^ I Michaelis-Menten equation still is an excellent approximation other than the two limits. 

o 
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The Michaelis-Menten (MM) mechanism [ij is widely used to understand the catalytic 
activities of various enzymes. According to this mechanism, a substrate S binds reversibly 
with an enzyme E to form a complex ES. ES then undergoes unimolecular decomposition 
to form a product P, and E is regenerated for the next cycle. 

E+S v^ ES -^ E° + P, E° -^ E (1) 

k-i 

The rate of product formation v on substrate concentration [S] can be characterized by the 
MM equation [l| 

where Vmax = k2[E]'Y is the maximum generation velocity, [E]t = [E] + [ES] is the total 
enzyme concentration, and Km = {k-i + kp)/ki is the Michaelis constant. Although the 
MM mechanism and equation have been found for almost a hundred years, they are still 
widely accepted and remain pillars of enzymology. 

Even if the classic MM equation achieves considerable success, there are still many in- 
triguing problems about the equation waiting to be answered. Particularly, the recent single- 
molecule fluorescence studies 2|, y, y, M, m| found that catalytic rates of many enzymes are 
fluctuating with time due to conformational fluctuations. A natural question hence is why 
MM equation works well despite the broad distributions and dynamic fluctuations of single- 
molecule enzymatic rates. Recently Xie et al. tried to address this issue from view of points 
of single-molecule experiment l7l| and theory jq] . In addition that the reciprocal of the first 
moment of f{t), {t)~^ = f/[-E']T follows MM equation well at any substrate concentration, 
the most remarkable discovery of their experiment is that the waiting time distributions 
f{t) exhibit highly stretched multiexponential decays at high substrate concentration and 
monoexponential decays at low substrate concentration [7[ . Xie et al. |8| attributed the non- 
exponential decay of f{t) to dynamic disorder of the rate constants of the reactions in Eq. (^ 
caused by transitions among different enzyme conformations. They theoretically proved that 
the classic MM equation still holds at the single molecular level when the transition rates 
among the ES conformations are slower than the catalytic rate k2 (the quasi-static condi- 
tion), even if f{t) is no longer monoexponential decays at high substrate concentrations. 
Therefore one of following issues is whether we can still derive the MM equation beyond the 



quasi-static disorder. Xie et al. 



indeed attempted to give an answer about it. But their 



effort ended in the two-state model for the algebraically complex in the multistate model. In 
this work, we propose that the classic MM equation holds under broader disorder conditions. 
Different from the discrete state model of Xie et al, a continuum diffusion- reaction model 
is used 9|. 

The conformational probability density for each enzyme state, Pi{x,t), in Eq. IQ can 
be obtained by three coupled diffusion- reaction equations with the potential [Vi(x)] and the 
reaction terms [A;j(x)] 

— PE(a;,t) = [Ce - hsix)] Pe + A;-i(x)Pes 

— PEs(a;, t) = [Ces - h{x)] Pes + kis{x)PE (3) 

d 

— PEo{x,t) = CeoPeo + k2{x)PES 

where 

Ci = Di— exp[-/3Vi(x)] — exp[/3Vi(x)], (4) 

and I=E, ES or E°, and k^{x) = A;_i(x) + k2{x) and kis{x) = ki{x)[S] are defined for 
convenience. The diffusion coefficient Di determines the rate of the conformational transition 
on the state i. The initial conditions are Pes{x,0) = 0, Peo{x,0) = 0, and Pe{x,0) is 
the thermal equilibrium distribution with the potential Ve{x). In single molecule turnover 
experiment, the observation is the probability density of the waiting time for an enzymatic 
reaction f{t), which is defined 

f{t) = I hix)PEsix)dx. (5) 

We first study the solutions of Eq. in two limiting cases: the slow reaction and the 
nondiffusion limits. In addition that the coupled diffusion-reaction equations have exact 
analytical solutions under these limits, this study would be useful in understanding the 
general solutions of Eq. Q . Particularly, we will show that quasi-static condition proposed 
by Xie et al. is just one of case of the latter limit. 

The slow reaction limit In this limit the processes of reactions is very slowly compared to 
processes of the enzyme conformational diffusion. Therefore the thermal distributions are 
always maintained during the courses of reactions. The solution to the diffusion-reaction 
equations then can be written as 

p,ix,t) = prix)p,it). (6) 



where Pi^{x) oc exp[— /3Vi(x)], and (3 ^ = k^T, k^ is the Boltzmann's constant, and T is 
absolute temperature. Substituting them into Eq. (^Hj) and considering that 

APr(x) = 0, (7) 

after simple calculations we get 

fit) = PEsit) J k,{x)P^l{x)dx (8) 

1.15 i,2 

2A, L J ' 



where Aeq 



(^k +41)74 -/^il 4s.. '^' and Seq = -(^s. + 41)/2. Hence the 



-'cq -"^cq ' ' J-'cq J^^eq 



-Jcq' 



reciprocal of the mean waiting time is 

1 kl,JS] 



{t) [S]+M, 



eq 



(9) 



where Meq = (A^es^ + ^isc )/4cc- ^^ ^^^ ^^^ ^^^^ ^^ ^'^^^ rapid diffusion limit, Eq. is 
almost the same as the single molecule MM equation in the absence of dynamic disorder [2| 
except that the rate constants now are the mean values on their inner conformational coor- 
dinate. 

The nondiffusion limit {k^^ <^ PDi) In this limit the reactions in Eq. (^ proceed so 
rapidly at the initial values of the slow coordinate x that the distribution of x is not restored 
by diffusion in the course of reactions. Then the diffusion terms in the diffusion reaction 
equations are neglected or Di ^ 0. The following calculations are simple and we immediately 
have 

f{t) = /'p^^(a:)^^l^M^ |e[sW+AW]* _ g[i?(.)-A(.)]t| ^^ (^q) 

J 2A[x) 



and 



(t) [S]+M, 



^ndPj /,,N 



nd 



where 



Mnd = Sid / P^'^{x)k3{x)/ [ki{x)k2{x)] dx, 



where A{x) = [(hsix) + hix)f/A - kisix)k2{x)f^ and B{x) = -{h + fcis)/2. 

We note that the expressions of the waiting distribution and the mean waiting time in 
the latter hmit is very similar with the main conclusion [Eq. (31)] obtained by Xie et al. J8|. 
It is not unexpected because the quasi-static condition used by Xie et al. is included in our 
nondiffusion limit. But two new points are revealed in the present work. One is that, in 
addition to /c2, the other rates may also be allowed to fluctuating in time. The other and 
more interesting finding is that the unknown steady-state weight function w{k2) introduced 
in prior by Xie et al. has a direct physical interpretation. To better understand the similarity 
between our calculation with them, we rewrite Eq. (fTUl) bv viewing /c2 as variable instead 
of X, and make use of the experimental observations [if that both ki{x) and A;-i(x) are 
independent of the conformational coordinate (the "wide reaction window" limit termed by 
Sumi and Marcus (ll|), then we obtain the same expression as that Xie et al. solved by very 
complicated algebra operations. 

m = [^(k-)^-^ [^''^'''' - ^''"'''] dk2, (12) 

while the weight function ^(^2) is related to the initial equilibrium distribution as follows, 

wik,) = P^"^ [x-\k,)] ^^, (13) 

where x~^{k2) is the inverse function of k2{x). In order to demonstrate the usage of this 
"microscopic" interpretation, we fit our theory by assuming that the potential Ve has a 
harmonic form with spring constant k, i.e., 

^E(x) = -^exp(-|^) (14) 

V27ra 2cr^ 



where a"^ = k^T/k, and k2{x) = aexp{—bx). It might be the simplest model of Eq. ()12p. 
The values of the parameters and fitting results are showed in Fig. ^ We see that our 
calculation is satisfactory. Because Eq. (fT^ is almost the same with previous result IS], we 
are not ready to explain the general behavior of it afresh. In the following part, we will 
focus on the general solutions to the coupled diffusion-reaction equations. 



Firstly substituting 



Id 
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Pi{x,t) = gi{x)Qi{x,t) (15) 



into Eq. Q, where gi{x), 1= E, ES, and E° are related to the thermal equilibrium distribu- 
tions 



- ^-/3Vi(x)/2 / 



-m{^)d, 



.1/2 



X 



we transform the diffusion reaction equations into an adjoint form 

d 

— gE(a;,t) = -[i^E + hsix)]QE + k'_^ix)QEs 
at 



d_ 

dt 
d 



QEs{x,t) = -[Hes + hix)]QEs + Ksix)QE 



-^QEo{x,t) = -HeoQeo + k2ix)QES, 

where the new functions k'_^{x), k'ig{x) and ki^lx) are respectively defined by 

k'_i{x) = k'_^gEs/9Eix), 
k[s{x) = kisgE/9Esix), 
k'^(x) = k2gEs{x)/gEoix), 



(16) 



(17) 



(18) 



and the Hamiltonian operators are 



Hi = -D, 



dx"^ 



PDj 



(3 fdVi 
2 V rfx 



d^Vi 
dx"^ 



(19) 



respectively. We assume that the operators Hi have discrete eigenfunctions \n)i (the bound 
diffusion assumption), i.e., 



Hi\n)i = ei„|ra)i, 



n = 0,l, 



(20) 



then gi{x) are just the lowest order eigenfunctions |0)i in the coordinate representation with 
zero eigenvalues (e^p = 0). The reader is reminded that the diffusion information has been 
included in the eigenvalues, for instance, given the potentials Vi to be harmonic like Eq. (|14j) . 
then ei„ = nk(3Di. Defining Oi = s + Hi + ki{x), here i = IS and 3 respectively correspond 
to I=E and ES, and Oeo = s + Heo, the Laplace transform solution Qes{x, s) of Eq. (fT7|) 
with the initial conditions can be written as 

1 



1; ' 



QEsix,s) = O^gk 



IS 



nrio). 



Oe ^-1^ES^15 



(21) 



Although the above calculations are exact formally, we cannot say more the inverse 
operator Ogg. Therefore we employ the decoupled approximation |llLll2l| 

l^kl^~'\0)u{0\k„ (22) 

where fcj = i(0|/i;j|0)i. This is would be exact when the expectation value of the operator 
Eq. (j^ is computed in the state |0)i. Using the approximated unit operators in Eq. (jSJ 
repeatedly, we finally get the analytical form of /(s) as follows, 

.. ^ _ E(|fc3fcls|0)Efc|s.,/fck, , . 

^^'' s' [1 + 4s(^)] [1 + ^iis)] - E{0\hhs\OhEs{0\k^ihs\0)Es/kilkh., ^ ^ 

where 

a\is) = kl^-\{0\h (^s + H,y' \0), (24) 

oo 

= kl^s-' + kl^~'Y,i^ + e,J-'\i{Omn),\\ 

n=l 

and i = IS and 3 correspond to I=E and ES, respectively. We immediately see that the 
waiting time distribution f{t) has a multiexponential behavior, because the denominator of 
Eq. ()23j) is a higher order (> 2) polynomial. For instance, if we truncate the aj(s) to nth, 
f{t) should be a sum of 2(n+l) exponential decay functions. A remarkable conclusion is 
that, even if the waiting time distribution function has very complicated multiexponential 
decay behavior, the reciprocal of the first moment of distribution, (t) = —df{s)/ds \s=o still 
has a simple MM-like expression, 

1 IC[S] 



(t) [S]+M' 
where 



(25) 



M = ki,JJ^, 

^ = kL (k'kl-^s{0\kik.,\0)ESE{0\k,ks\0)E /kLklS / J'k'Em,ks\0)E 



OO 



and 

^ = ( 1 + kl,^^-'Y.'EL |Es(0|A;3|n)Es|' ) 4., + 4.,"' E^e„' |E(0|A;i|n)E|' A;^ 

\ n=l / n=l 

Here we have separated the substrate concentration [S] from the rate k^c^jx). Under the two 
limiting cases discussed at the beginning, Eq. (j^ is approximated to be 



xll 



a](s) ^ kl^ s ^ (slow reaction limit), 

a]{s) ~ ki{x)s~^ (nondiffusion limit). (26) 



Substituting them into Eq ()23p and making the Laplace transformation, we obtain the same 
Eqs (jHl) and pO|) . The general solution hence well recovers the two limiting cases. Because 
the decoupling approximation Eq. (j^ has been proved to be a good approximation [ll|, |l2| , 



we conclude that classic MM equation still is a good approximation even in the presence of 
dynamic disorder with arbitrary characteristics. 

There are two main contributions in the present work. Firstly we recover the waiting 
time distribution f{t) obtained by Xie et al. in quasi-static condition, and given a 
microscopic interpretation of the weight function used by them. But compared to their 
complicated algebra calculation and a continuum approximation involved, our approach 
is very simple and direct. We must point out that the current calculations except fitting 
to the experiment are independent of specific conformational dynamics. Second, we get 
general waiting time distribution Eq. ()23|) with arbitrary dynamic disorder, and prove 
that the reciprocal of its first moment still follows the classic MM equation. Although 
this conclusion is based on decoupling approximation, it still is positive because this 
approximation has been proved to work well in various systems. Moreover, it is beyond the 
quasi-static disorder condition. While the discrete chemical reaction scheme of Xie et al is 
hard to achieve because of mathematic difficult. We believe that Eq. ()23|) would be more 
useful than Eq. (fTUj) when experiments are performed on various enzyme molecule under a 
broad range of environmental condition. 
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FIG. 1: Waiting time distribution vs. substrate concentration. The dotted and dashed Hnes and 
the experiment data are from Ref. 0. The substrate concentrations are 10 fiM (the cross), 20 
/iM (the circle), 50 /iM (the time) and 100 /uM (the square), respectively. The parameters used in 
Eq. (O are A;i = 5 X lO^M^^s^^ fc_i = 18300 s"\ a = 904 s'^, b = 5.0 nm and a = 0.1 nm [ij. 
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